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Uncertainties in cloud feedback remain in GCMs 


Local contribution to intermodel spread 
in cloud feedback: AR4 


Local contribution to intermodel spread 
in cloud feedback: AR5 



# * 


' O’ 




Regression of Local Cloud Feedback to Global-mean Cloud Feedback (W-rrv^K 1 ) 

Most of intermodel spread arises from low stratocumulus/cumululs regions 



• Low subtropical clouds still uncertain. 
Soden and Vecchi (2011 • Large contribution from equatorial Pacific. 


Soden and Vecchi (2011): 

• Low cloud cover is responsible for ~3/4 of the difference in global-mean net 


cloud feedback among AR4 models, with the largest contributions associated 
with low-level subtropical marine cloud systems; 





CAM5, CAM5 (IPHOC; CLUBB), and AM3 

(CLUBB, CLUBB+) 



Model state is updated 
after each process 
(“sequential-split”) 
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Model state is updated 
after each process 
(“sequential-split”) 
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Advance 12 prognostic equations 
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Use PDF to close higher-order 
moments, buoyancy terms 


Select PDF from given family 
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Differences between IPHOC and CLUBB used in 

GCMs? 

CLUBB (Cloud Layers Unified by Binormals; Golaz et al. 2002); 
IPHOC (Intermediately Prognostic Higher-order turbulence 
Closure; Cheng and Xu 2008) 
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Global Distribution of Annual Mean Low Cloud Fraction - 

- IPHOC 


a)CAM5 mean= 41.7 rms= 1 5.45 corr= 0.68 
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Global Distribution of Annual-Mean SW Cloud-radiative 

Forcing - IPHOC 


a)CAM5 mean= -50.6 rms=16.18 corr= 0.86 


b)CAM5-IP mean= -47.9 rms=14.48 corr= 0.85 




Global mean from CAM5-IP is the closest to CERES 
slightly weaker negative forcing from low clouds 






CAM5, CAM5-CLUBB (tuned) cloud fraction 


and SW cloud radiative forcing 
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camel ubb_F2000 (yrs 1-5) 

TOA SW cloud forcing meon= -48.92 
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GFDL AM3, AM3-CLUB and tuned versions 


SW Cloud radiative forcing differences from CERES 


(a) AM3-MG: tune (NA in Table 2), 
v=l,no Enhance_Accr 


(bl) AM3-CLUBB:no tune,v=l,no Enhance_Accr 



60E 1Z0£ ISO 120W 60# 

Ion - (0..360.) lot - {~9Q..aO.) r(Oba, Mod) - 0.853592 

Ma-d — Obs =■» -5.37706 rms* = 13.32Q5 


GOE 1 2C£ 180 120# 60# 

Ion - (Q, F 36Q r ) l at - <-90„90.) r(Obg. Mad) - 0.765392 
Mad - Obs - -16.SB65 rm 5fl - 24 5505 


Parameter 

Tuned 

Original 

Cl 

1.0 

2.5 

C4 

1.0 

5.2 

C5 

0 

0.3 

C6 

0.5 

4.0 


(b2) AM3-CLUBB: tune only 
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(b3) AM3-CLUBB: variable v only 
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(b4) AM3-CLUBB: tune + variable v (b5)AM3-CLUBB:tune+variable v+Enhancc. 
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GDFL AM3 united parameterization, CLUBB+ 
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Tuned parameter tests in CAM5-CLUBB (Guo etal. 

2015) 


Table 1 . Tunable 
Para meter 

Parameters of CLU BB and ZM 

Description 

Default Value 

Investigated Range 

Cl 

Constant associated with w 2 dissipation 

2.5 

1.25-5 

C2rt 

Constant associated with q t 2 dissipation 

1.0 

0.5-2 

C6rt 

Low skewness of Newtonian damping of water flux 

4.0 

3.0-8 .0 

C6rtb 

High skewness of Newtonian damping of water flux 

6.0 

3.0-8 .0 

C7 

Low skewness of bucyancy damping of water flux 

0.8 

0.25-1.0 

C7b 

High skewness of buoyancy damping of water flux 

0.65 

0.25-1.0 

C8 

Constant associated with Newtonian damping of w 3 

5.0 

1. 5-6.0 

CI1 

Low skewness of buoyancy damping of wj_ 

0.8 

0.0-1 .0 

Cl 1b 

High skewness of buoyancy damping of iv 3 _ 

0.65 

0.0-1 .0 

C14 

Constant of Newtonian damping of u 2 and v 2 

1.0 

1 .0-2.0 

v (nu) 

Background coefficient of eddy diffusion 

20.0 

10.0-40.0 

fi [beta] 

Constant related to skewness of (?i and qt 

1 .75 

0.0-3.0 

y (gamma_coef) 

Constant of the width of PDF in w -coordinate (er^) 

0.52 

0.1 -0.6 

fi (mu) 

Parcel entrainment rate (1/m) 

0.001 

0.5-2.0 X 10 " 5 

COJnd 

ZM precipitation efficiency overland 

0.0059 

0.003-0.09 

C0_ocn 

ZM precipitation efficiency over ocean 

0.045 

0.005-0.09 

dmpdz 

Entrainment rate of ZM 

-10 -5 

-0.2 to -2 X 10 3 

tau 

CAPE consumption time scale (s) 

3600 s 

1800-10,800 




Sensitivity to Tuning parameter tests in CAM5- 
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Summary and conclusions 

• The higher-order turbulence closure approach offers a promising 
approach to subgrid-scale variability. 

• The low-level clouds are improved in different GCM simulations and 
the biass in SW cloud radiative forcing are reduced. 

• The potential for realistic simulation of cloud processes is great with 
the higher-order turbulence closure approach, for example, coupling 
with cloud microphysics, and unified low and deep convection 
parameterization. 

• Sensitivity to parameters are especially strong for skewness-related 
parameters. A better constraint is needed from global observations. 



